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Abstract 

We compare the measurement of the gravitational mass of 38 high-redshift 
galaxy clusters observed by Chandra using Modified Newtonian Dynamics 
(MOND) and standard Newtonian gravity. Our analysis confirms earlier 
findings that MOND cannot explain the difference between the baryonic 
mass and the total mass inferred from the assumption of hydrostatic equi- 
librium. We also find that the baryon fraction at r25oo using MOND is 
consistent with the Wilkinson Microwave Anisotropy Probe (WMAP) value 
of Q B /Q M . 
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1. Introduction 



Since the introduction by Milgrom in 1983 (Milgrom 1983), MOND 



had success in explaining galaxy rotation curves using only the mass-to- 
light ratio as a free parameter, and was able to predict the applicability of 
the Tully-Fisher relation to low surface brightness galaxies before dynamic 
information on them was available (Scarpa 2006). These successes are not 



surprising, given the fact that MOND was created as a phenomenological 
theory in order to eliminate the need for dark matter in galaxies. 

While there have been numerous applications of MOND to individual 
galaxies, it is also important to study MOND on those scales where dark 
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matter is believed to dominate, in particular that of galaxy clusters. There 
have been far fewer studies of MOND on this scale, with previous work on 
the subject including Sanders] (fl999|) , Aguirre et al. (20011, Pointecouteau 



k Silk| p)05| ) and |Angus et al.| Q2008[ ). These results indicate that MOND 
does not eliminate the need for dark matter in galaxy clusters. In this paper 
we use MOND to calculate gravitational and baryonic masses for a sample 
of 38 galaxy clusters. The motivation of this work is the need to confirm the 
previous results by using a larger sample of massive clusters at high redshift 
(z=0.14-0.89). 

We begin with a brief overview of the data and the models used, found 
in Section [2} Section [3] describes the calculation of the gas mass using X-ray 
observation, the derivation of the MOND acceleration, and the calculation 
of total masses from hydrostatic equilibrium. We then present our results 
in Section [4] and our conclusions in Section [5] The cosmological parameters 
h = 0.7, f2M=0.3 and f2A=0.7 are used throughout this work. 



2. Chandra X-ray Data And Data Modeling 



We analyze Chandra X-ray data from 38 clusters in the redshift range 
z=0.14 — 0.89 (Table [T]). The Chandra observations and the data model- 



ing with the isothermal /3 model are presented in Bonamente et al. (2006); 



LaRoque et al. (20061; Bonamente et al. (20081. Here those aspects of the 



data modeling and analysis that are relevant to the present investigation. 
The electron density model is based on the isothermal spherical /?-model 



(Cavaliere & Fusco-Femiano 1976 1978), which has the form 



n e (r) = n eQ 1 + 



-3/3/2 



(1) 



where n e o is the central electron number density, r c is a core radius, and (5 
is a power-law index. The radial profile of the X-ray surface brightness is 
obtained via integration along the line of sight, and results in the following 
analytical expression: 



,r0 



1 + 



2\ (i-6/fl)/a 



(2) 



Best-fit model parameters and confidence intervals for all model param- 
eters are obtained using a Markov chain Monte Carlo (MCMC) method 



described in detail by Bonamente et al. (2004) and LaRoque et al. (2006). 



For each cluster, the Markov chain constrains the parameters S x o, /3, 
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T e , and the chemical abundance (see LaRoque et al. 2006 for best-fit val- 
ues). We use the cosmological parameters h = 0.7, fiM=0.3 and U\=0.7 to 



calculate each cluster's angular diameter distance Da (e.g., Carroll et al. 



1992). 



All of our calculations are done out to a maxium radius of r25oo> the 
radius at which the cluster mass density is 2500 times the critical density, 
which is also approximately the radius out to which our Chandra data are 
sensitive without any extrapolation of the models. In LaRoque et al. (20061 
we have compared masses obtained using the simple isothermal (3 model 
with those obtained using a more complex non-isothermal f3 model, and 
shown that the two methods yield the same ratio of baryonic to total mass. 
The non-isothermal model included a double (3 model distribution for the 
density (Equation 7 in LaRoque et al. 2006), and an arbitrary temperature 
profile that was constrained by assuming that the plasma is in hydrostatic 
equilibrium with an NFW potential (Equation 8 in LaRoque et al. 2006). 
We therefore expect no significant bias in the study of MOND masses using 
this simple model, which has the advantage of analytical expressions for the 
observables and the masses. 



3. Calculation of gas, Newtonian and MOND Masses 



For the isothermal [3 model, the gas mass enclosed within a given cluster 



radius is given by (e.g. LaRoque et al. 2006 ) 



M gas {r) 



A 



r/D A 



1 + k) 



9 2 \ "3/3/2 



8 2 d8, 



(3) 



where A = 4tt fi e n eo m p D A , [i e is the mean electron weight (calculated from 
the Chandra data, with typical value of [i e = 1.17), and Da is the cluster 
distance. The central electron density n eo is calculated from the X-ray 



surface brightness (LaRoque et al. 2006 Birkinshaw et al. 1991) as 



n,, 



S x Jtt(1 + z) 4 (/W Me )r(3/?) 
A eH D A 7r^(3P - l/2)9 c 



1/2 



(4) 



in which hh is the mean hydrogen weight (typical value fin = 1-4), and A e # 
is the emissivity of the plasma. The gas mass in Equation ^ is measured 
directly from observables and does not depend on the law of gravity in any 
way. 
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For the total mass, Mtotah we solve the hydrostatic equilibrium equation 
using Newtonian gravitation and the ideal gas law, to obtain the Newtonian 
mass: 



M ma i(r) 



kr 2 



Gn e (r)ntotm p 



dn e (r) ,dT e (r) 
T e (r)— hn e (r) 



dr 



dr 



(5) 



in which fi to t is the mean molecular weight (typical value fitot = 0.61). Owing 
to our use of an isothermal model, the second term goes to zero and using 
Equation ([I]) it becomes 



Mtatalir) 



3/3kT e 



G[im v r~ + r 2 



(6) 



Starting with the most basic form of MOND (Pointecouteau & Silk 



2005; Scarpa 20061, we derive the gravitational mass based on the MOND 



law of gravity. The MOND acceleration (a) is related to the Newtonian 
acceleration (ajv) by the relation 



GM{r) 



a N = a - fi( — ). 

a Q 



(7) 



where a is a critical acceleration that is meant to be a new physical constant 
and represents the acceleration below which Newtonian gravity is no longer 



valid (Milgrom 1983 Scarpa 2006); the function ^ is any interpolating 
function that provides a transition between the two regimes. We take the 
critical acceleration to be ao =1.17e _8 cm/s 2 as used in Angus et al. (2008); 



this value is derived from galaxy rotation curve fits Scarpa ( 2006 ) . In Section 



[4] we also discuss the results obtained using a variable value of the critical 
acceleration. 

The interpolating function fi(x) is any function which yields the appro- 
priate asymptotic behavior, 



H(x) 



1 for x 1 
x for xCl 



The most used interpolating function (Scarpa 2006; Pointecouteau & Silk 



2005), and the one which we use here, is 



x- 



(8) 
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Using Equation Q in Equation (J7| and solving for a, gives us the relation 

--'11 r~n\ 



'M2 + 2\/ 1+ S-° 



JV 



Equation d9|) is the MOND acceleration as a function of the Newtonian 



acceleration and the critical acceleration parameter (see also Pointecouteau 



& Silk 2005). We use the hydrostatic equilibrium equation according to 



MOND, 



1 dP n 



a(r), 



(10) 



p g (r) dr 

{p g {r) is the gas density and P g its pressure) in order to determine the 
MOND equation for total mass analogous to Equation ([6| for Newtonian 
gravity. Using Equations (fot) and (10) with = GM(r)/r 2 leads to the 



elimination of the Newtonian acceleration a at: 



GM{r) 



-1 dP 

p(r) dr 



(11) 



al + 



-1 dP 

p(r) dr 



For the j3 model profile of density we obtain: 



3f3KT 



M M ond{t) 



prrip r 2 c + r 2 



(12) 



„ / 3/3KT r 
a£ + 



V pm p r£ 



rj. + r l 



This equation can be further simplified when one realizes that the numera- 
tor is the total mass according to Newtonian dynamics, Equation (|6]). This 
simplification gives us a final equation for the total mass based on the hy- 
drostatic equilibrium condition and using MOND gravity. 



M M ond(t) 



M Newt0 n{r) 



i 



( a pm p (rj + r 2 ) 
^ 3f3KT r 



(13) 



This equation displays the appropriate asymptotic behavior, namely that 
for a critical acceleration of zero we recover the Newtonian Mtotai- 
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4. Results 



4-1. Comparison of Newtonian and MOND gravitational masses 

Having found an equation for the MOND total mass (Mmond) we cal- 
culate M gas , Mtotal an d Mmond for each of the 38 galaxy clusters in our 
sample out to a radius of r25oo- In each case M gas was calculated using 
Equation ([3]), Mtotal by means of Equation ([6]), and Mmond using Equa- 
tion (13). For each cluster, the MOND gravitational mass is lower than the 
the Newtonian gravitational mass by just ~ 20 — 30%, and both are signif- 
icantly larger than the gas mass M gas (see Table [T]) . The radial profiles of 
the masses are shown for all clusters in Appendix A. 

In order to calculate the baryonic mass we approximate the stellar con- 
tribution by M* = 0.15M 9as , as is done in Pointecouteau & Silk (2005). 
We therefore estimate the total baryonic mass as Mb ar yon = M gas + M*. 
We find that the average f gas {r250o) for MOND is 14.2±0.2%, and that 
the average fbaryon(r 2 5o o) is 16.4±0.2% (Table |2~1). T hese results agree with 
the results obtained by Pointecouteau & Silk (2005), which report a value 
of /fearyon=20.2±2.0% at a larger radius of riooo- We find that, even with 
the inclusion of M*, the baryon mass is still significantly lower than the 
total MOND mass for every cluster, confirming the fact that MOND is 
incapable of eliminating the need for dark matter in clusters. The mean 
f gas and f baryon in Table [2] are calculated as the weighted mean of all mea- 
surement, in which the weight is the standard error of each measurement. 
We also calculate the un-weighted means and their root-mean-square er- 
rors as / 9as =0.109± 0.003 (Newton) and / 9as =0.139 ± 0.004 (MOND), and 
fbary On =0- 125 ±0.004 (Newton) and f bary on =0.1 60 ± 0.005 (MOND). 

We show the average distributions of fbaryon an d M to t a i/MMOND as func- 
tion of physical radius in Figure [TJ The distribution of fbaryon shows that 
the dynamical mass is always significantly higher than the baryonic mass. 
The distribution of M tota i/ Mmond shows that the difference between the 
Newtonian and MOND masses are largest near the cluster core and at the 
largest radii, where the gravitational accelerations are lowest and therefore 
the MOND correction largest. Typical values of r25oo are between 0.5 and 1 
Mpc, therefore we do not extrapolate these plots beyond 1 Mpc. The Chan- 



dra data presented in this paper were fit with a 100 kpc core cut (LaRoque 



et al. 2006), and therefore we do not extrapolate our models below this 



radius either. 

The cluster baryon fraction calculated via the standard Newtonian dy- 



namics in this paper {{fbaryon) = 0.1305) and by other authors (e.g., Mc- 



Carthy et al. 2007 for a review) are systematically lower than the WMAP 
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Table 1: Cluster Parameters And Masses For The Isothermal /3-Model 



Cluster 



?"2500 

(arcsec) 



201.6ffJ 



(Gpc) 



M gas 
(1O 13 M ) 



2.63^ 



Mtotal 

(1O 14 M ) 



9 ,rW 

z - lo -0.14 
+0.36 
-0.36 
-0.5 
-0.42 

4.82±g;l 5 

2.49l°i! 

4.96+ - 8 
3.32 



MMond 

(1O 14 M ) 



+0.14 
-0.12 
+0.32 



Abelll413 

Abelll689 

Abelll835 

Abelll914 

Abelll995 

Abell2111 

Abell2163 

Abell2204 

Abell2218 

Abell2259 

Abell2261 

Abell267 

Abell370 

Abell586 

Abell611 

Abell665 

Abell68 

Abell697 

Abell773 

CLJ0016+1609 

CLJ1226+3332 

MACSJ0647.7+7015 
MACSJ0744.8+3927 
MACSJ1149.5+2223 
MACSJ1311.0-0310. 
MACSJ1423.8+2404 
MACSJ2129.4-0741 
MACSJ2214.9-1359 
MACSJ2228.5+2036 
MS0451.6-0305 . . 
MS1054.5-0321 . . 
MS1137.5+6625. 
MS1358.4+6245 . 
MS2053. 7-0449 . . 
RXJ1347.5-1145. 
RXJ1716.4+6708 
RXJ2129.7+0005 
ZW3146 



0.143 
0.18 
0.25 
0.17 
0.32 
0.23 
0.2 
0.15 
0.18 
0.16 
0.22 
0.23 
0.38 
0.17 
0.29 
0.18 
0.26 
0.28 
0.22 
0.54 
0.89 
0.58 
0.69 
0.54 
0.49 
0.55 
0.57 
0.482 
0.41 
0.55 
0.83 
0.78 
0.33 
0.58 
0.45 
0.81 
0.24 
0.29 



+5.1 
-5.5 



217.1 
171.1±| 4 
226.7+ 4 - 6 
133.4+ 4 -^ 



+9.8 
-8.9 



140.5 
206.7lt 9 



256.11 
190^ 
172 
145 
131.3 
97.95 



-13 
-11 
!+5.7 
-5.4 

172. 1+^jj 
0+6:7 

•°-6.3 
-8.5 
-7.5 
+4.1 
-4.1 

181.6±f 

no.3+l;| 

160 - 5 -33 
hlO' 
-9.3 

133.ll| 9 
148.9+|^ 
79.85±1 



153' 



66.02 



+7.4 
-6.5 



91.94±|;| 
58.88+g 
70.64+l:| 



+7.6 

-6.6 



73.92 
65.93+j- 1 
72.42+ 4 -? 
91.311^ 



83.61 
82.18 



-4.1 
-3.8 
+3.7 
-3.5 



33.74+7+ 
41.63+|:| 



113.4 



-5.8 



54.2 



5.2 
+5.1 

4.6 



122.3+J-jj 

45.01+t 4 

+5.3 
4.9 



128.5 



131 



?+2.6 
•°-2.5 



0.52 
0.63 
0.81 
0.6 
0.96 
0.76 
0.68 
0.54 
0.63 
0.57 
0.73 
0.76 
1.07 
0.6 
0.9 
0.63 
0.83 
0.88 
0.73 
1.31 
1.6 
1.36 
1.47 
1.31 
1.25 
1.32 
1.35 
1.23 
1.12 
1.32 
1.57 
1.54 
0.98 
1.36 
1.19 
1.56 
0.78 
0.9 



+0.17 
-0.17 
f0.23 



5.1 

5 - 81 -0 2 

4.87+0.12 
3.51 



0.1 

+0.14 
0.14 



2 19+ - 27 

z - ly -0.23 

8-08+Jf 
474+0.27 

3.02 



0.24 
+0.13 
0.12 



1 Q9+0.14 

1 - 82 -0.14 
o nQ+0.2 

2- 24±g : ? 8 
2-78±8:l 1 

o 97 +0.13 
z - z ' -0.11 
9 07+O.H 
z -°' -0.11 

2 63+ - 1 

z - oo -0.095 
o «c;+0.34 
°- DO -0.31 

^•^-0.26 

2.74+ - 17 



-0.15 

4.38+0.29 
3.89" 



-0.29 
1+0.51 
-0.46 

4.911^ 
3.08^:1 
3.09±°| 4 

9 -,9+0.25 
z - ±z -0.21 

2 26+ - 1 

z - zu -0.096 
q 97+O.3 
°- z ' -0.26 

O4+0.32 
°- y ^-0.31 

2^70+0.23 
•' 8 -0.21 

4.8l°,f 
0.9051°:^ 

1 26+ - 15 

i - ZD -0.14 
hO.19 
-0.17 
-0.13 



2.53; 
0.933 



-0.12 
123 +o.i9 
2.57 H 



'-0.16 
7+0.16 
-0.14 

4.44±°:^ 



4.95 
5.55 



0.64 
+0.31 
0.27 
+0.27 



l-7 y -0.24 

2 ^6+ - 36 
Z - OD -0.31 

2 03+ - 42 
Z - UO -0.33 
9 70+O.37 
z - ' °-0.33 

9 A O+0.34 

z - 4 °-0.28 
9 -1 q+0.22 
z, - lo -0.2 
9+0.14 
z -0.12 



4 S2+ - 91 
q 4fi +0.4 

2 59+ - 3 

qq + 0.39 

°-°°-0.37 

5.971}-? 
2 26+ - 42 

z - zu -0.35 
9 q -1+0.4 
Z - O1 -0.34 
9 1 7+0.74 
z - 1( -0.53 

1 Q4+ - 2 

9 09+O.57 
z - 8z -0.45 
q oc-+0.68 
°- oo -0.59 
h0.32 
-0.27 
-0.53 
-0.46 

61 2+ - 47 
u - Diz -0.36 

1 02+ - 26 

1 - uz -0.2 

q 1 q + 0.51 

°- lo -0.41 
-I 99+O.38 
± - zz -0.29 

8 - 08 ia69 

1 QO+0-45 
1 - oy -0.33 

2.08l°:i 
3.62l°: 22 



2.051 
3.76" 



1.51 
3.97 



-0.33 
4.59+0 1 : 4 ! 

q or;+0.27 
°-°°-0.22 
q m+0.46 
°- yl -0.44 
1 O/I+0.48 
i - 84 -0.37 
4 47+0.22 
-0.21 
q qr;+0.72 
°- yo -0.57 
9 trq+0.27 
z.oo_ 24 

-I 9 q+0.22 
1 - zo -0.19 
1 q+0-31 
± - y -0.26 
1 4fi +0.35 
i - 4D -0.27 
9 iq+0-33 
z - ly -0.29 
1 o+0.29 
llO -0.23 
1 co+0.19 

1 -^ 1 -0.099 

3-48±g;i 
9 70+0.36 

z -'°-0.32 
1 OQ+0-26 
i.yo_ 23 

9 o+0. 36 
z -°-0.34 

4.81H 

5.3111' 2 

1.89l°:l 

1.86+of 
1 7-1+0.66 

i -0.47 
1 c 4 +0.17 
1 - O4 -0.16 
9 qrr+0.53 
z -°°-0.41 
q 9/1+0.63 
°- Z4 -0.54 
1 c-7+0.28 
i - ' -0.23 
q 9+0.49 
°- z -0.43 

0.454l° : 44 

0.802l°- 23 
1-0:45 

-0.36 
h0.33 
-0.24 
7 1 o+0.73 
'• i8 -0.65 
1 14+0-41 

1 c;+0.22 
1 -°-0.19 
9 07+O.I9 
z -°'-0.18 



2.47: 
0.921 



Table 2: Newtonian And MOND Mass Fractions 



Cluster 



/ {gas} 

Newton 



MOND 



n 19 o+o.oo4i 
u.izo_ 0045 

ins+ - 0045 
u.iuo_ 0041 

+0.0047 

-0.005 

+0.0034 

-0.0037 

h0.0051 

-0.0046 

n nsso+o.0078 

n -.a 7 +0.0026 
u - i4 ' -0.0025 

n ncK7+ - 0094 

U - Uy57 -0.0099 



-0.008 



f{baryon} 

Newton 



MOND 



1-0.0048 



Q 2Q1 +0.0092 
u.zui_ 0098 

n 1 /IS+0.0078 
U ' 148 -0.007 

U.l^D_g g 082 

o i46+ - 0058 

n mo+o.0084 
u.±uo_ 0074 

n ioy+0.017 
u - 10 ' -0.016 

n ons+0-0047 
0.z08_ 00046 

n i qq+o.016 

U - iJ8 -0.016 
n -, 07+O.OO8I 
U- 1 - 3 ' -0.0082 
n i 7 +0.016 
' -0.015 

n isq+0.017 

n -1 77+O.O23 
u - ±l '-0.022 

U.14D_ 00096 
n 14 r+0.014 

n 17 o+0.012 
U-J-'^-O-Oll 
n 91 c+0.0083 
v.z,±o_ 00g5 

n 1 9 +0.016 
U - 1Z -0.014 

o-i86ig:8g 
o.i63i8;8l? 
o-isl^i 1 
o.o93H8:8S 

u.iuo_ 0012 

U - i8 ' -0.019 
n 101 +0-015 
u - ly± -0.014 

0.14218:81 

0.16818:81} 

n ifil +0.019 

u.ioi_ 0018 

n 1 4 +0.015 
U - i4 -0.013 

20S+ - 018 

0.17218:81^ 
0.21918:8^ 
o.i8l8:8l 
o.n8l8:811 

116+ 0023 

U.±±O_ 02 
n 1/19 +0.0082 

1 24+ - 022 

n 107+0.017 

u - ly '-0.016 
n 1 7S+0-0076 
u - 1 ' o -0.0073 



Abelll413 

Abelll689 

Abelll835 

Abelll914 

Abelll995 

Abell2111 

Abell2163 

Abell2204 

Abell2218 

Abell2259 

Abell2261 

Abell267 

Abell370 

Abell586 

Abell611 

Abell665 

Abell68 

Abell697 

Abell773 

CLJ0016+1609 

CLJ1226+3332 

MACSJ0647.7+7015 
MACSJ0744.8+3927 
MACSJ1149.5+2223 
MACSJ1311.0-0310. 
MACSJ1423.8+2404 
MACSJ2129.4-0741 
MACSJ2214.9-1359 
MACSJ2228.5+2036 
MS0451.6-0305 . 
MS1054.5-0321 . 
MS1137.5+6625 
MS1358.4+6245 
MS2053.7-0449 . 
RXJ1347.5-1145 
RXJ1716.4+6708 
RXJ2129.7+0005 
ZW3146 



0.175 



0.105 

o.ior 

0.0741 



0.0909 



0.131: 

0.0843 
0.127 
0.106 
0.1311 
0.0746 
0.0821: 
0.136' 



Weighted Mean 



h0.0042 
-0.0044 
n ln9 +0.0069 
u - luz -0.0069 
n 1 1 c+0.0084 
U.1±0_q 0082 

u - ill -o.on 

n 1 +0.0054 
u - 1 -0.0052 

0916+ 0068 

n 1 1 1 +0.0058 
U.- 1 11 -0.0057 
0.0036 
0.0038 
h0.0089 
-0.0083 
0.0066 
-0.0063 
+0.0058 
-0.0055 
+0.0067 
-0.0063 
-0.015 
-0.014 
0.0091 
0.0084 
0.012 
0.012 
1 o, 4 +0.0077 
u.io^±_ 0078 

0.0979l8'81 7 

n 1 1 fi +0.0'064 
Uillo -0 .0061 

u - ilD -0.011 
n i n9 +0.0088 
u.±uz_ 0084 

n iqc+0.0087 
u.ioo_ 0087 

fl 197+0.0088 
U - 1Z ' -0.0085 
n i 4 r+0.015 
U.14O_ 012 

o.i23i8:81t 
0.080918:88^ 
o.o76ii8:81? 

-0.0056 
L -0.0055 

0.0883l8on 
hO.o'osa 

-0.0077 

n 1 90+0.0044 
u.izj_ 0043 



o.ir 
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Figure 1: Average distribution of fbaryon and M to tai/MMOND for the 38 
clusters; the dashed lines represent the root-mean-square uncertainty of the 
38 measurements at each radius 



value of Ob/Om=0.176±0.02 (Bennett et al. , 20031, indicating the presence 



of undetected baryons in clusters. Within MOND, the results of this paper 
{{fbaryon) = 0.164) and those of Pointecouteau & Silk (20051 are in statis- 



tical agreement with the WMAP value of Qb/^m- We therefore speculate 
that MOND can in principle reconcile the cluster baryon fraction with the 
cosmic value of Qb/^m, without the need for additional undetected baryons 
in clusters. 



4-2. Results for a variable value of the MOND critical acceleration 

We also consider the possibility of ao being a free parameter of the 
MOND theory, and use our data to determine it. To this end, we begin 



with Equation (13) for the MOND total mass, and calculate the value of the 
critical acceleration at each radius such that Mmond{i~) = Mb aryon {r), 



■aryon 



(r) = M M ond{t) 



M Newt0 n{r) 



(14) 



\^3PKT 



I a {im p r 2 + r 2 



We then solve equation ( 14 ) for a and find that 



o (r) 



2.I3KT 



lim p + r z 



-^Newton 
baryon 



1, 



(15) 



9 



where M^ ewton is given by equation Q. The radial distribution of a is 
shown for all clusters in Appendix B. 

In Figure [2] we show the value of the critical acceleration calculated at 
^2500 for all 38 clusters, with a mean value of o = 1.49±0.09 x 10~ 7 cm/s 2 . 
The main result is that the value of the "free" a at ^2500 differs from that 
measured from galaxy rotation curves (e.g., Scarpa 20061 by one order of 
magnitude. Moerover, we find statistically significant differences in the mea- 
surement of a Q at different radii for a given cluster (see Appendix B), and 
for different clusters at the same radius of r25oo (Figure [2j . 




cluster 



Figure 2: Calculation of a (^250o)- The horizontal line is the mean of a Q 
and its root-mean-square error, calculated as the un-weighted mean of the 
38 measurements. The weighted mean of the 38 measurements is a =1.17± 
0.02 x 10~ 7 cm/s 2 . 



We also show the average values of the best-fit critical acceleration a 
for the entire sample, as function of radius, in Figure [3j At all radii, the 
critical acceleration required to explain the hydrostatic MOND mass is at 
least one order of magnitude larger than the canonical value obtained from 
galaxy rotation curves (e.g., Scarpa 2006). 



5. Conclusions 

In this paper we have measured gravitational and gas masses at ^2500 
for a sample of 38 high-redshift galaxy clusters using standard Newtonian 
gravity and Modified Newtonian Dynamics. This is the largest sample of 
galaxy clusters to which the MOND theory of gravity is applied to date. 
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Distance(kpc) 

Figure 3: Average distribution of a Q for the 38 clusters; the dashed lines 
represent the root-mean-square uncertainty of the 38 measurements at each 
radius 



We initially used a fixed value of the MOND critical acceleration a = 
1.17e~ 8 cm/s 2 and measured an average (fbaryon) = 0.164; this result con- 
firms that MOND does not eliminate the need for dark matter in galaxy 
clusters. Our work therefore confirms the previous reports by Angus et al. 



Q2008H and |Pointecouteau k Silk] ( |2005[ ) . Further evidence that the MOND 
theory of gravity is unable to explain X-ray observations of galaxy clusters 
is found by a fit of our X-ray data to a MOND model with varying a Q . The 
value of a required to achieve a baryon fraction of unity varies both with 
cluster radius, and between clusters; moreover, the mean value of this free 
a is one order of magnitude larger than that required by galaxy dynamics. 
We therefore conclude that X-ray observations of the hot cluster plasma do 
require strong presence of dark matter even when MOND is used. 

This analysis also finds that the fbaryon at r25oo calculated using MOND 
is in statistical agreement with the WMAP value of 0{,/Omj similar to the 
results of Pointecouteau & Silk ( 2005 1 . 
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Appendix A 

In this appendix we provide the radial profiles of the gas, Newtonian and 
MOND masses for all clusters, calculated between 100 kpc and r25oo- 

Figure 4: Distribution of M gas ,M t otah and Mmond as a function of radius(arcsec) for 
all 38 clusters. The symbols are used as follows: O f° r M to tai, x for Mmond, and A for 

M gas 
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(19) Abell773 (20) CLJ0016+1609 (21) CLJ1226+3332 




(22) MACSJ0647.7+7015 (23) MACSJ0744.8+3927 (24) MACSJ1149. 5+2223 




(25) MACSJ1311. 0-0310 (26) MACSJ1423.8+2404 (27) MACSJ2129.4-0741 
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Appendix B 



In this Appendix we provide the radial distribution of the best-fit values 
of the critical acceleration an, for all clusters. The plot extend between 
100 kpc and r25oo- 

Figure 6: The critical acceleration a as a function of radius, under the condition 
Mmond = Mbaryon for all 38 clusters 
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